Lens design based on instantaneous focal function 
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The formula for the lens is derived based on the information of instantaneous focal function. Focal 
function is an important tool in designing lenses with extended depth of focus (EDoF) because this 
allows EDoF lens designers to try out various mathematical curves using computers to optimize 
their design. Once an optimal focal function information is obtained, the corresponding physical 
EDoF lens can be fabricated using the lens equation formulated in this presentation. 



I. INTRODUCTION 

Optical imaging system with large depth of field 

a)OF) is required to produce sharp images (author?) 
. In photography, the DOF is the portion of a scene 
that appears sharp in the image, for example, the region 
denoted by A in Fig. [TJ Ordinarily, a lens focuses paral- 
lel rays of light at one distance known as the focal point, 
as illustrated in Fig. [21 Therefore, not all points within 
the DOF can be claimed as focused per se. However, 
due to a gradual decrease in the sharpness of the im- 
age from the focused spot, the amount of blurring within 
the DOF is imperceptible to human eyes under normal 
viewing conditions. As such, in particularly for films and 
photography, the image region can be subdivided into 
two, where one lies within and the other lies external to 
the DOF. In the photograph of Fig. [U the region A, 
wherein the image appears sharp and focused, is said to 
lie within the DOF; whereas the region B, in which re- 
gion the image is blurred, is said to lie external to the 
DOF. 




by either decreasing the size of lens aperture or by in- 
creasing the shutter speed, or through tweaking of the 
both. These methods, however, drastically reduces the 
amount of light passing through the lens and require ex- 
tra lighting. For developing still images, the required ex- 
tra lighting can be accommodated by the use of a fiash. 
For motion pictures, however, this approach proves to be 
inadequate, as typical video involves about thirty frames 
of images per second, [l^ In addition to this difficulty, the 
smaller lens aperture increases diffraction and this places 
a practical limit on the extent to which the DOF of an 
optical imaging system can be enhanced by the afore- 
mentioned methods. That being said, can the DOF of 
an imaging system be increased (A) without sacrificing 
the intensity of light passing through the lens and (B) 
without increasing the diffraction? The answer to this 
question is yes and this approach involves some sort of 
digital filtering. 
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Figure 1: Illustration of depth of field. Region A is within 
the depth of field (DOF), whereas the region B is outside of 
DOF. The image is blurred drastically in region B. 



One focal point, |3 



Figure 2: An ideal ordinary lens is characterized by a single 
focal point. 



Alternatively, but equivalently, the DOF in an imag- 
ing system is defined as the distance in the object space 
in which objects are considered to be in focus. The dis- 
tance over which objects appear sharp can be increased 
by extending the DOF of an imaging system. Tradition- 
ally, the DOF of an imaging system can be increased 
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The digital filtering method requires a scheme for the 
image reconstruction algorithm based on the principles 
of wave optics(author?) 0, S, 0, H, S, 0| • The image re- 
construction code is often hardcoded in the accompany- 
ing digital processing unit (DPU) and, for this reason, 
the imaging method based on digital filtering is coined 
as the software assisted imaging technology or SAIT for 
short. lis] The SAIT solution for an imaging system is a 
promising technology in that it has potential to increase 
the DOF while encompassing altogether the processes of 
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which (A) drastically reduce the amount of light passing 
through the lens and (B) those problems associated with 
increased diffraction due to the reduced diameter of the 
lens aperture. 

The software assisted imaging technology is not as per- 
fect as it sounds and it has some of its own problems to 
offer. Among them, one directly relates to the image 
processing speed. In the traditional imaging solutions 
based on lenses, which I refer to as the analog technology, 
complex arrangement of lenses function to focus image 
at the focal plane from wherein the image is developed. 
As such systems based on analog technology by-pass the 
digital filtering stage altogether, images are developed in- 
stantaneously in analog imaging systems. However, due 
to the aforementioned problems and limitations of the 
traditional method for increasing the DOF, the analog 
technology provides only a limited solution when con- 
cerning the image quality. The systems based on SAIT 
principle face problems that are exactly the opposite in 
nature from that of the analog systems. In principle, the 
reconstructed image by a way of image reconstruction al- 
gorithm can be made to resemble the original image to 
any level or degree of resemblance, provided the image 
processing speed is of no concern. But, such methodol- 
ogy would Hmit systems based on SAIT to still imagery 
applications and the video sector of the market must be 
discarded, which is a bad idea for business. Naturally, 
for systems based on SAIT principle, a trade off must 
be made between the image quality and the image pro- 
cessing speed. 

Among the early pioneers to successfully commercial- 
ize imaging system based on SAIT are Cathay and 
Dowski, who did much of their work at the University 
of Colorado (author?) 0, 0, 0]. In the modern liter- 
ature, their work is cited as "wavefront coding." The 
SAIT solution based on wavefront coding has been trade 
marked by CDM Optics, Inc., and it is known as the 
Wavefront Coding"^". The imaging solution based on 
wavefront coding basically involves two stages: the in- 
put and the output stages. The input stage involves the 
optical element and this represents the hardware contri- 
bution side of the SAIT solution. The output stage in- 
volves the DPU, wherein the image reconstruction code 
base is hardcoded, and this represents the software con- 
tribution side of the SAIT solution. The optical element 
in SAIT solution is distinguished from the input stage 
of traditional imaging system, which is just complex se- 
ries of ordinary lenses, in that it produces many focal 
points along the optical axis instead of just one at the 
focal plane. Such optical element is referred to as lens 
having extended depth of focus (or EDoF for short) and 
this is illustrated in Fig. [31 

Both the quality of reconstructed image and the image 
reconstruction speed are critically important in SAIT. In 
principle, the ordinary lens, such as the one illustrated 
in Fig. [21 can just as well serve as the input stage for the 
SAIT imaging systems. However, this must be done at 
the cost of overly complicated algorithm routines for the 



software side of the system. The length of image recon- 
struction code directly relates to the number of transis- 
tors in a DPU. As a simple rule of a thumb, more tran- 
sistors there are in a DPU, more energy it requires to 
operate. And, more lines of coding for the image recon- 
struction algorithm implies the slower processing speed 
for the reconstructed image output. The full-HD qual- 
ity of a video involves sixty image frames every second. 
This implies, the application of SAIT system to full-HD 
video processing would require the image reconstruction 
time span of 16 ms or less. The demand for very fast 
processing speed and low power consumption make ordi- 
nary lenses inadequate for the input stage of the SAIT 
system, which leave open for an alternate solution for the 
optical element to be used as the input stage. 
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Figure 3: Lens with extended depth of focus has many focal 
points, /i, /2, /a, and so on. An ideal EDoF lens has infinitely 
many focal points and all light rays are confined within the 
cylindrical tube very small diameter w and very long length 
I. 



The image reconstruction code base can be optimized 
for image processing speed and quality if the EDoF lens 
is used as the optical element for the input stage of SAIT 
system. The EDoF lens is characterized by parameters 
I and w. The parameter / represents the depth of focal 
points along the optical axis and the parameter w repre- 
sents the width of bundle containing light rays associated 
with each focal points, as illustrated in Fig. [31 The ideal 
EDoF lens has the parameter I of which is infinite in 
length and the parameter w of which is infinitely thin. 
For the realistic EDoF lenses, however, the parameters 
I and w are typically on the order of microns. 

The idea of EDoF lens as an optical element which fo- 
cuses Hght into longitudinally directed fine along the op- 
tical axis was first proposed by Golub, et. al. (author?) 
Q. The idea was later adopted by others and it has found 
applications in various imaging systems, such as micro- 
scopes, cameras^and lithography to Hst a few, (author?) 
0, [3, Si S [1; [Hi- Alexander and Lukyanov have 
recently proposed a conceptual scheme for the EDoF 
lens(author?) jTl|. Their scheme for EDoF lens con- 
sists of zones that are axially symmetric about the opti- 
cal axis and this is illustrated in Fig. [H The idea behind 
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their concept is as follows. The light ray crossing each 
zone gets focused to a unique spot on the optical axis 
and this spot is within in the EDoF lens parameter I. 
If Pi is the function which describes the focal length for 
the ith concentric zone in Fig. [H the lens aperture, in 
principle, can be tailored to behave like EDoF lens by 
tweaking Borrowing their terminology, the (3i is re- 
ferred to as the instantaneous focal function. The idea 
behind EDoF lens is to make I as large and w as small as 
possible. By experimenting with different functions for 
Pi, the EDoF lens parameters I and w can be engineered 
to the acceptable range for the imaging systems based on 
SAIT principle. This is exactly just what Alexander and 
Lukyanov did, and their focal function is summarized in 
Fig. m In the figure, the vertical Hues represent discon- 
tinuities and each of the twelve zones has been indicated 
appropriately by a number. For each of the zones in Fig. 
m the incident parallel rays are focused at different points 
on the optical axis within I; and, the focusing is described 
by the instantaneous focal function, P = f3i, summarized 
in Fig. m 
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Figure 4: Schematic of conceptual EDoF lens proposed by 
Alexander and Lukyanov. The lens aperture has an axial 
symmetry about the optical axis. 



input from a lens aperture satisfying the focal function 
described in Fig. [H the physical EDoF lens to be de- 
signed had a constraint of satisfying the same focal func- 
tion characteristics. To end the story, the physical EDoF 
lens with such characteristics for the focusing behavior 
was found(author?) The obtained physical profile 
of lens was entered into CODE V@[l3| to generate the 
needed PSF information for the lens aperture. [l^ This 
information was in turn used by Alexander to test for the 
performance of his algorithm. 



Instantaneous focal length (P) vs. lens radius (y) 
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Figure 5: Focal function /3 proposed by Alexander and 
Lukyanov (author?) [H] ]. 



This work concerns the result of my role in the project, 
which was to design a physical EDoF lens from the focal 
function characteristics. As such, this work concerns the 
hardware contribution side of the optical imaging system 
based on SAIT. 



The imaging system based on SAIT is a package 
solution in which both the hardware (optical element 
or EDoF lens) and the software (image reconstruction 
code base) contributions must be optimized. Alexander 
and Lukyanov experimented with various mathematical 
curves for the focal function P in an attempt to make 
the EDoF lens parameters I and w as ideal as possible, 
but did not attempt to provide any solutions concern- 
ing the physical shape or the profile for their conceptual 
EDoF lens. To test the image reconstruction code for 
the processing speed and the quality of generated image 
output, the information on the point spread function (or 
PSF for short) of the optical element input stage is re- 
quired. Since no information on the physical profile of 
their conceptual EDoF lens was available, the test for 
their image reconstruction algorithm had to be deferred. 
It was my job to design a physical EDoF lens. Since 
Alexander's image reconstruction code was based on the 



II. THEORY 

A. Axial symmetry 

By definition, an axial symmetry is a symmetry about 
a given given axis. The object has an axial symmetry if 
its appearance is unchanged with the rotation about some 
axis. Illustrated in Fig. [6] is a schematic of conceptual 
lens, which shows an axial symmetry about the optical 
axis, i.e., the x axis. Such a lens can be dissected through 
the origin with the xy plane as shown in Fig. [6l The 
curve traced on the xy plane, which is a set of points 
on the surface of lens, can be revolved about the optical 
axis for the three dimensional shape of the lens. The 
design of axially symmetric lens, therefore, simplifies to 
the problem of finding the set of points on the surface of 
lens of which gets traced on the xyplane. 
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Lens cun'ature on xy plane 



Figure 6: Schematic of lens symmetric about optical axis. 

The physical law of which governs the bending of light 
is the Snell's law. I shall apply Snell's principle to derive 
the equation for the cross-sectional profile of the lens. 
The terminology, "cross-sectional profile of the lens," im- 
plies the lens curvature on the xy plane, which is illus- 
trated in Fig. [6l 



B. Derivation from Snell's law 

When a ray of light passes across media of different 
refractive indices, its path is governed by the Snell's law. 



710 sm ( 



ng smfy. 



(1) 



as illustrated in Fig. [71 Here, = (cu) and ng = 
ng (ui) are frequency dependent refractive indices with uj 
denoting the angular frequency of the light. The param- 
eters (j) and 6 represent the angle of incidence and angle 
of refraction, respectively. 

If N denotes the normal vector to the local point y — j 
on the curve x = h{y) , then it can be shown 

li-N X (-ei)ll = ll-NII ll-eill sin = TV sin 

and the expression for sin becomes 

. , IINxeil 



N 



iV^IINII, 



(2) 



where ei is the unit basis for the x axis. 

Similarly, the expression for sin 9 may be obtained by 
considering vectors A, B, and C of Fig. [7l The vectors 
A, B, and C satisfy the relation, 

A + B^C. (3) 

In explicit form, vectors A and B are defined as 

A = -762, B = (/3-a)ei, (4) 

where 62 is the unit basis for the y axis. With Eqs. ([3]) 
and the vector C becomes 



The vector cross product N x C is given by 

N X C = (/3 - a) N X ei - 7N X 62 

and its magnitude becomes 

UN X C|| = - a) N X ei - 7N X e2|| 

= NC sine, (6) 

where N = ||N|| and C = ||C|| . Utilizing Eq. jS]), C may 
be expressed as 



C= (C-C)'/' = ^{i3-af+r 

and the Eq. ([6]) is solved for sin 9 to yield 

|(/3 - a) N X ei - 7N x 62 



1/2 



smf = 



TV 



1/2 



(7) 



Refracted ray 



X^h(y) 



'angent line 




Figure 7: Illustration of Snell's law. The scalar quantity (3 
denotes the focal length. 

Insertion of Eqs. ^ and ([7]) into the Snell's law of Eq. 
(HI) gives 



n0 
ng 



(/3 — a) N X ei - 7N x 62 1 



Nxeill (/3-a)^+72 



1/2- 



(8) 



{j3-a)ei -762. 



(5) 



By definition, the normal vector N satisfies the relation, 

g{x,y) ^x-h{y), 

where g {x, y) is a function whose gradient gives N, 

„ dg dg dh 
N = V5 = — ei -I- —62 = ei - —62. 
ox oy ay 



5 



Because N is the normal vector at the location 
(x = a, ?/ = 7) , I write 



9h 

N = ei - — 
dy 



62- 



(9) 



y=7 



The following vector cross products are valid, 
N X ei = ei X ei 

N X 62 = ei X 62 



dh 
dy 
dh 
dy 



62 X ei, 



62 X 62, 



where Eq. ([9]) was used to replace N. Since ei x ei 
62 X 62 = 0, the previous relations reduce to 



9h 

N X ei = — 

dy 



63, Nxe2 = e3, 



(10) 



■y=7 



where 63 is the unit basis for the z axis of which satisfies 
the relation, 

ei X 62 = 63, 62 X ei = -63. 

Insertion of Eq. (fTO|) into Eq. ^ gives 



ne 



dy 



y=i 



■ 7 



dh 
dy 



1/2 ' 



which expression can be rearranged to yield 

dh 
dy 



13 — a 



1/2 ' 



where a and 7 are constants of which are depicted in Fig. 

El 

For Alexander and Lukyanov's optical element, the in- 
stantaneous focal function (3 = (3 (y) in Eq. ifTTj) is as 
defined in Fig. [H The 7 for the y axis is not anything 
special, of course. Any y belonging to the domain of h 
satisfies the Eq. ifTTj) . The generalization of Eq. (fTTj) for 
all y belonging to the domain of h is done by making the 
following replacements: 



a ^ X, 7 



dh 
dy 



V=7 



dh dx 
dy dy' 



With these replacements, Eq. ifTTj) gets re-expressed in 
form as 



dx 

dy 



y 



(3-x- 



"0 



{(3 - xY + y2 



1/2- 



(12) 



How is the instantaneous focal function, /?, restricted? 
The 13 in Eq. I|12p is restricted so that the expression for 



dx/dy does not blow up. Equation (fT2)) is well defined if 
and only if the denominator satisfies the condition. 



P-x-'^Up-xf + yA^'%Q. 
ne L J 



(13) 



Contrarily, but equivalently, the previous statement can 
be reworded as follows. Equation lfT2|) is ill defined if and 
only if the denominator satisfies the condition. 



f3-x-^ {13-xY + y'' 
ne L 



1/2 



0. 



(14) 



For reasons to follow, I shall proceed with the latter. To 
solve for /3, I shall first rearrange Eq. lfT4|) as 



13 -x='^ \{l3-xf + y'' 
ne L 



1/2 



(15) 



Squaring of both sides give 



{(3-x) ^ ^{(3-x) + -^y 



This expression can be rearranged to become 



= 32 y 



or 



{13 - xY 



(11) And, solving for {(3 — x) ,\ obtain 



2 

2 

2y ■ 



(3-x = ± 



n<t,y 



(16) 



where the ± came from the action of taking the square 
root on both sides, of course. Now, one of the signs in 
Eq. (fTBl) can be eliminated by comparing with Eq. (flSl) . 
This is the reason why I proceeded with Eq. ([M]) instead 
of Eq. ifTS]) . The and ne in Eq. ifTS]) are both real re- 
fractive indices, which cannot be negative numbers. The 
instantaneous focal function, (3, and the lens thickness, x, 



{P-xY + y^ 
1/2 



must be non- 



must be real, which implies 

negative else {[3 — xY + y^ becomes an imaginary 

term. As real refractive indices cannot be negative num- 
bers, the term {(3 — x) is also a non-negative real in Eq. 

(flS)) . provided n^ > 0, ne > 0, and {(3 — xY + y^ > 0, 

of course. Therefore, the right hand side of Eq. fT6|) 
must be positive; and, this gives 



l3^x + 



T'4>y 
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Now, this is precisely the condition for (3 which makes Eq. 
(fT2l) ill defined. Equivalently, then Eq. lfT2|) becomes well 
behaved for (3 satisfying the condition given by 



(17) 



Equation ifTTl defines the restriction for the instanta- 
neous focal function, f3. 

What can be concluded of the restriction so defined in 
Eq. lfT7|) for the instantaneous focal function? To answer 
this, recall that terms such as (3, x, y, n^, and ng are all 
real values. And, there are no restrictions on and ne 
to speak of which of the two must be bigger or smaller in 
value. Interesting per se, the choice of > ng results in 
the statement, 



> ng, 



(18) 



where the i denotes the imaginary symbol and the term 

— in Eq. (flTl) has been modified to • 

But, this condition defined in Eq. fT8|) is always satisfied, 
as /? is a real function. Therefore, it is concluded that 
Eq. (fT2|) is well behaved everywhere for > ng. 

III. RESULT 
A. Lens surface equation 

The profile of axially symmetric lens about its optical 
axis is obtained by solving the initial-value differential 
equation, Eq. (fT2| . 



dx 
dy 



y 



13-x- 



(/3 - xf + y^ 



1/2 ' 



x{yo) = xo, 



where x (yo) = xq is the initial condition to be specified 
and the instantaneous focal function /3 satisfies the con- 
strain defined in Eq. lfT7|) . Without loss of generality, one 
may choose x {y = yo = 0) = for the initial condition 
and the lens profile satisfies the differential equation. 



dx 

dy 



y 



P-x-^Y-P-xf + y^ 
n4,y 



1/2 ' 



(19) 



a;(0) = 0, (3^x + 



The quantities and ng are the two refractive indices in 
which one represents the lens and the other representing 
the surrounding medium. Which of the two refractive 
indices corresponds to the lens depends on the configura- 
tion of the problem, as demonstrated in the proceeding 
sections. 



B. Alexander and Lukyanov lens 

1. Instantaneous focal function 

The instantaneous focal function proposed by Alexan- 
der and Lukyanov has been discussed previously in Fig. 
[5l The instantaneous focal function for each of the twelve 
zones can be curve fitted and represented by a quadratic 
polynomial of the form given by 



P = Pi ^ ayf + byi + c, 



< 



(20) 



where the subscript i of {Pi, yi, yi,mm, yi.max) denotes the 
i th concentric zone. The coefficients a, b, and c, and the 
range for y, which defines the width for each of the axially 
symmetric concentric zones, are summarized in Table [H 
Since the instantaneous focal function, P, and the lens 
radius, y, have units of length measured in meters [m] , 
the coefficient a must have a unit of [m~^] , c a unit 
of [m] , and b must be a unit-less scalar. To reduce the 
width of the table, couple columns were represented in 
millimeter units, [mml . 



Table I: Domain yt and coefficients (a, b, c) of (3i — ayf + byi + c 
for Fig. \5\ 



yi.min, J/i.max [mm] 


a[l/m] 


b 


c [mm] 


0.0, 0.19182692 


-313.07 


0.0235 


3.5137034 


0.19519231, 0.27259615 


534.53 


-0.2472 


3.527877626 


0.27596154, 0.33317308 


-309.02 


0.0818 


3.5088062 


0.33653846, 0.38701923 


536.05 


-0.3275 


3.5493232 


0.39038462, 0.43413462 


-306.12 


0.1182 


3.502912672 


0.4375, 0.47451923 


539.03 


-0.3891 


3.569538239 


0.47788462, 0.51153846 


-303.68 


0.1463 


3.496845208 


0.51490385, 0.54855769 


542.21 


-0.4417 


3.589312176 


0.55192308, 0.58221154 


-301.27 


0.1695 


3.49080193 


0.58557692, 0.6125 


545.96 


-0.4895 


3.609151039 


0.61586538, 0.64278846 


-298.81 


0.1893 


3.484870978 


0.64615385, 0.67307692 


179.08 


-0.0474 


3.469596542 



2. Lens profile 

Equation fT9|) was solved by the Runge-Kutta routine 
coded in FORTRAN 90(author?) [H]. For the com- 
putation, refractive indices, n^ and ng, were chosen as 
follows, 

710 — 1.5311, ng = 1.0. 

In this configuration, n^ denotes the refractive index for 
the lens and ng denotes the refractive index for air.[l9| 
The resulting cross-sectional profile of the lens curvature 
is shown in Fig. [HI The profile in Fig. [8| was revolved 
about the optical axis, which is the x axis in the figure, 
to generate the three dimensional profile of the physical 
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lens. This result is shown in Fig. O where the optical axis 
is located at {Xy = 200, Xz = 200) . The scaling factor of 
A — 3.365385 x 10^^ was introduced for graphing purpose 
only. 

In spite of the non constant (3 for the instantaneous fo- 
cal function (see Fig. [5]), the resulting EDoF lens shown 
in Figs. [8] and [9] seems to resemble the parabolic curve, 
which configuration is known to have only one focal point, 
i.e., (3 = constant. (author?) Is the result portrayed 
in Fig. [8] (or Fig. [9]) correct? To give a qualitative an- 
swer to this, I shall recall the EDoF lens parameter I, 
which was previously illustrated in Fig. [31 

Cross-sectional profile of lens curvature 



^ -4C-05 




«e= 1 



where /3niax and /3min denote, respectively, the maximum 
and the minimum values in the profile of instantaneous 
focal function. In the focal function profile of Alexander 
and Lukyanov, Fig. \S[ /?max ~ 3.5188 mm and /?min ~ 
3.4831mm. Plugging the information into Eq. l(2T|) . this 
roughly gives / « 36 um. The resemblance of the EDoF 
lens to the paraboHc curve can be attributed to the small 
value for I. Considering that Alexander and Lukyanov's 
lens has a maximum radius of y « 0.67 mm, which can 
be identified from Table [H the lens diameter comes out 
to be about d = 2y fa 1.34 mm. This implies, the lens 
diameter is larger than the EDoF lens parameter, I, by a 
factor of thirty seven. Under such circumstance, the lens 
could be perceived as having a single focal point from 
the perspective of human eye. In spite of the existence 
of number of very closely spaced focal points within the 
length of / along the optical axis, the human eyes do not 
have sufficient resolving power to distinguish those focal 
points. Even less so, the human eyes cannot distinguish 
the actual cross-sectional profile of the lens of which is 
only sHghtly perturbed from the cross-sectional profile of 
the paraboHc lens. To prove that this is indeed the case, 
I shall use the derived lens formula, Eq. lfT9|) , to generate 
various parabolic lenses. 



0.0002 0.0003 0.0004 
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Validation of the result 



Figure 8: The cross-sectional profile of lens curvature cor- 
responding to the instantaneous focal function proposed by How does one know that Eq. 1191 generates the cor- 
Alexander and Lukyanov, Fig. \5\ rect profile for the lens? The easiest way to settle this 

dilemma is to actually apply Eq. (fTg)) to the well known 
types, i.e., the parabolic lenses. 

Tln'ee dimensional piofik of leus 
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Figure 9: Three dimensional profile of lens satisfying the focal 
function /3 proposed by Alexander and Lukyanov, Fig. [5] The 
y and z axes have been multiplied by A = 3.365385 x 10~^. 

In explicit form, the EDoF lens parameter / is defined 

as 

^=l!/3max-/3min|l, j ^""^ I!' (21) 



1. Simple parabolic lens 

The parabolic curves are known to merge parallel 
rays of incidence light to a unique focal point called 
a focus (author?) [lj|. As a consequence of this, the 
parabolic curves fiatten in the curvature with the focus 
positioned at distances further from the vertex. Such 
property of parabolic curves make them ideal for testing 
and validating the lens formula defined in Eq. (fT9l) . The 
constant focus of /3 = 1 m, /3 = 5 m, and P — 10 m were 
considered to generate curves using Eq. (fT9| : and, the 
result is summarized in Fig. [Tol As expected, the gen- 
erated curves are that of parabolic curves in which each 
curves corresponds to focal points = 1 m, /3 = 5 m, and 
P = 10 m. The curve corresponding to /3 = 10 m is more 
fiat in curvature than the ones corresponding to /? = 1 m 
or /3 = 5 m, as expected. The result corresponding to 
(3 — Im was revolved about the optical axis to illustrate 
the three dimensional profile of the parabolic lens. This 
is shown in Fig. [iTJ Again, in the figure, the optical axis 
is at {Xy = 200, Xz — 200) , where the scaling factor A is 
A = 3.365385 x 10"^ 
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Cross-sectional profile of lens curvature 
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Figure 10: The cross-sectional profile of lens curvature for 
/3 = 1 m, /3 = 5 m, and P = lOm. 



tion are, respectively, /Smax = 2.5 m and P^iin = 0.01m, 
which can be verified from Table HIl Using the formula 
for I defined in Eq. l(2T|) . this gives I w 2.49m. Therefore, 
in this test configuration, I is much larger than d, i.e., 
d w 5.4 X 10~^Z, which is just the opposite situation from 
that of Alexander and Lukyanov. That being said, the 
Eq. (flOl) was plotted for the curve and the result is shown 
in Fig. [131 As expected, the resulting cross-sectional pro- 
file for the lens does not resemble simple parabolic lens. 
However, the curve profile for each of the twelve zones 
in Fig. [13] represents the portion of a paraboHc curve 
corresponding to Pi illustrated in Fig. [121 By superim- 
posing the two graphs. Figs. [12] and [131 the boundaries 
for each zones can be identified by kinks in Fig. [13] The 
cross-sectional profile of the lens curvature illustrated in 
Fig. [T3l was revolved about the optical axis for the three 
dimensional profile of the lens. This result is shown in 

Fig. m 



Three dimensional profile of lens Instantaneous focal length (P) vs. lens radius (y) 




V, [m] **venjc![J lines are diseontinuilies** 



Figure 11: Three dimensional profile of a lens with a constant Figure 12: Focal function /3 with characteristic of a step func- 
instantaneous focal function, f) — Im. The y and z axes have tion. 
been multiplied by A = 3.365385 x 10"^. 



2. Concentric parabolic lens 

Now I consider a slightly more complicated profile for 
the instantaneous focal function, (3. To show that, in- 
deed, the likeHness of Alexander and Lukyanov's EDoF 
lens to the parabolic lens is attributed to the small value 
for the EDoF lens parameter I, I shall modify only the 
/3 portion of Alexander and Lukyanov's profile for the in- 
stantaneous focal function, while leaving the size of lens 
diameter unmodified. The modified focal function for 
this test is shown in Fig. [12] And, the coefficients a, 6, 
and c, and the range of y corresponding to the (3 = [ii oi 
Eq. l(20l) for each of the twelve zones is summarized in 
Table [ll] Besides the increased focal length for each of 
twelve zones, the focal point is unique within each zones 
in this test configuration. In Alexander and Lukyanov's 
focal profile, / was much smaller than the lens diameter 
d, i.e., d ~ 371. The /3max and /3min for the test configura- 



Table II: Domain yt and coefficients (a, b, c) of (3i — ay^ 4- 
byi + c for Fig. [12] 



i/i,min, yi.max [mm] 


a[l/m| 


6 


c[m| 


0.0, 0.19182692 


0.0 


0.0 


0.25 


0.19519231, 0.27259615 


0.0 


0.0 


0.05 


0.27596154, 0.33317308 


0.0 


0.0 


0.50 


0.33653846, 0.38701923 


0.0 


0.0 


0.10 


0.39038462, 0.43413462 


0.0 


0.0 


1.00 


0.4375, 0.47451923 


0.0 


0.0 


0.05 


0.47788462, 0.51153846 


0.0 


0.0 


2.50 


0.51490385, 0.54855769 


0.0 


0.0 


0.10 


0.55192308, 0.58221154 


0.0 


0.0 


0.50 


0.58557692, 0.6125 


0.0 


0.0 


0.10 


0.61586538, 0.64278846 


0.0 


0.0 


0.01 


0.64615385, 0.67307692 


0.0 


0.0 


1.00 
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Cross-sectional profile of lens curvature 
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Figure 13: The cross-sectional profile of lens curvature corre- 
sponding to the focal function defined in Fig. [121 



Three dimensional profile of lens 

(- 
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Figure 14: Three dimensional profile of lens satisfying the 
focal function defined in Fig. 1121 The y and z axes have been 
multiplied by A = 3.365385 x 10"^ 

Basing on these results, it can be concluded that for 
the case where the lens diameter d is much larger than 
the EDoF lens parameter I, the profile of the lens re- 



sembles closely the profile of simple parabolic lens. Here, 
the word "simple" has been used to denote the parabolic 
curve with single focus. In the opposite situation, where 
the lens diameter d is much smaller than the EDoF lens 
parameter I, the profile of the lens no longer resembles 
the simple parabolic lens. Instead, in this latter case, the 
shape for the lens resembles superimposed, multiple num- 
ber of parabolic lenses of different degrees of curvature, 
as illustrated in Fig. [141 



IV. CONCLUDING REMARKS 

The image processing speed and the quality of pro- 
cessed images are both of critical importance in software 
assisted imaging technology. Such a requirement calls 
for the optimization of image reconstruction code based 
on the principles of wave optics. The coding side of the 
SAIT system can be optimized if EDoF lens is used for 
the input stage. 

In this presentation, the formula for the EDoF lens has 
been derived based on the knowledge of instantaneous fo- 
cal function, (3. The (3 information is an important tool in 
the design of EDoF lens, as this allows optical engineer 
to try out various mathematical curves using computers 
for optimization. With the knowledge of (3, this can be 
achieved without having to actually make EDoF lens 
prototypes, thereby saving time and the cost. Once the 
optimal solution for the instantaneous focal function is 
obtained, the physical EDoF lens can be manufactured 
based on the lens formula presented in this work. 
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